Viability of primordial black holes as short period gamma-ray bursts 
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It has been proposed that the short period gamma-ray bursts, which occur at a rate of ~ 10yr -1 , 
may be evaporating primordial black holes (PBHs). Calculations of the present PBH evaporation 



rate have traditionally assumed that the PBH mass function varies as Mj 



-5/2 



This mass function 



only arises if the density perturbations from which the PBHs form have a scale invariant power 
spectrum. It is now known that for a scale invariant power spectrum, normalised to COBE on 
large scales, the PBH density is completely negligible, so that this mass function is cosmologically 
irrelevant. For non-scale-invariant power spectra, if all PBHs which form at given epoch have a 
fixed mass then the PBH mass function is sharply peaked around that mass, whilst if the PBH mass 
depends on the size of the density perturbation from which it forms, as is expected when critical 

5/2 

phenomena are taken into account, then the PBH mass function will be far broader than M BH . 
In this paper we calculate the present day PBH evaporation rate, using constraints from the diffuse 
gamma-ray background, for both of these mass functions. If the PBH mass function has significant 
finite width, as recent numerical simulations suggest, then it is not possible to produce a present 
day PBH evaporation rate comparable with the observed short period gamma-ray burst rate. This 
could also have implications for other attempts to detect evaporating PBHs. 



PACS numbers: 98.80Cq, 98.70 Rz 



I. INTRODUCTION 

The physical origin of gamma-ray bursts (GRBs) is one 
of the outstanding problems in astrophysics Q . The ma- 
jority of gamma-ray burst are now known to have cosmo- 
logical origin, however Cline and co-authors have stud- 
ied the short (duration < 200ms) GRB population [Q, 
and suggest that they may be due to the evaporation 
of primordial black holes (PBHs) located in the galac- 
tic halo P]. These short events, which make up roughly 
2% of the total population, have simple time histories 
and hard spectra, relative to the longer duration GRBs, 
and are consistent with a Euclidean source distribution, 
suggesting a local origin 0. 

Primordial black holes (PBHs) are black holes which 
form in the early universe There are a number of 
possible mechanisms for their formation, including the 
collapse of cosmic string loops || and the collision of bub- 
bles formed at phase transitions The most natural 
formation mechanism is the collapse of large inflationary 
density perturbations Jt|— [Tc|] . PBHs evaporate via the 
emission of Hawking radiation jll] and PBHs with mass 
Mbh = M* ~5x 10 14 g would be evaporating today. In 
the standard picture of PBH evaporation (l2| all parti- 
cles with rest mass less than the black hole temperature 
Tbh where 
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dM BH "(Mbh) 



df 
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For the standard model of particle physics a (Mbh) > 
7.8 x 10 26 g 3 s- 1 for M BH - 5 x 10 14 g. PBHs 

of this mass, which are evaporating today, would have 
Tbh > 20MeV in the final stages of their evaporation. 
Cline and collaborators argue that q;(A/bh) could then be 
significantly increased by the effective number of degrees 
of freedom due to the quark-gluon phase transition [B|, 
resulting in a burst of radiation with duration consistent 
with the observed duration of the short period GRBs. 

To calculate the PBH evaporation rate the PBH mass 
function and local density enhancement are needed. 
Cline and collaborators take the bound on the global 
number density of PBHs per logarithmic mass interval, 
J\f g , at M = M*, from the diffuse gamma-ray back- 
ground to be < 10 5 pc~ 3 . They combine this with 
a local density enhancement factor rj = p\/p g = 5 x 10 5 , 
where T and 'g' denote local and global values respec- 
tively, to obtain the local number density of PBHs per 
logarithmic mass interval at M = M*, Af\ = rjJ\f s ~ 
10 10 pc -3 . The local PBH evaporation rate is then given 
by 
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and Mbh is the PBH mass in grams, are emitted. The 
rate of mass loss therefore depends on the number of 
particle degrees of freedom, a(M): 



which, using eq. (|2[) and the fact that evaporating PBHs 
can only be detected within ~ lpc 0], becomes 
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The constraint on the global PBH number density No = 
d^BH/dlnMeH has traditionally been calculated [ |14|Jl3| 
by assuming that the PBH mass function is a power law 
drtBH~/dA/BH oc Af B ^ 2 in which case the initial PBH 
number density is given by |L9] 



v(Muq) = 9.5 x 1CT 5 @ using the COBE normalisa- 
tion |19| . Using the Press-Schechter formalism Kim and 
Lee found the initial mass function produced by a power 
law power spectrum |2Cf|: 



^PBH.O Pc 
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dn B H 
dM BH 



n 



where fipBH.o is the present day fraction of the critical en- 
ergy density, p c in PBHs. Using MacGibbon and Carr's 
1991 evaluation of the diffuse gamma-ray bound (which 
assumes the M -5 / 2 mass function) p4j: 
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ilpBH < 7.6(±2.6) x 10 _9 /i 



-1.95±0.15 



(6) 



where h is the Hubble parameters in units of 
lOOkms^Mpc -1 , gives 7V g = 4.2 x 10 3 pc~ 3 . 

The power law PBH mass function, dnBH/dA/BH oc 



A/, 



-5/2 
BH 



where pi and Mm are the energy density and horizon 
mass when the PBHs form, immediately after reheating 
finishes at the end of inflation and the universe becomes 
radiation dominated. The horizon provides a sharp lower 
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cut-off in the mass function at Mbh = 7s Aifjji ■ 

For n ^ 1 the exponential term cuts the mass function 
off sharply so that it is peaked around A/bh ~ Mm- Only 
if n = 1 does the mass function have the traditionally 
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, was derived by Carr in the 1970s |l5[ for scale- used M BH 7 form. However for a scale invariant power 



invariant density perturbations. In this paper we out- 
line the arguments which show that this mass func- 
tion is cosmologically irrelevant. We then calculate the 
constraint on the global PBH evaporation rate for the 
sharply peaked mass function which arises from assum- 
ing that all PBHs which form at a given epoch have the 
same mass, and also for the far broader mass function 
found by recent numerical studies ]l6|] , due to near crit- 
ical gravitational collapse jL7|. Finally we review the 
calculation of the local density enhancement factor. 



II. PBH MASS FUNCTION 

Early studies by Carr [pl|, which assumed that all 
PBHs which form at a given epoch have the same mass 
(Mbh ~ 7s^ 2 Mh where Mh is the horizon mass at 
that epoch and 7 S parameterises the equation of state: 
p — 7s/o), found that an extended PBH mass function is 
only possible if the primordial power spectrum is scale- 
invariant (equal power on all scales). For Gaussian dis- 
tributed fluctuations the probability distribution of the 
smoothed density field p(S(Mn)) is given by 



spectrum normalised to COBE on large scales the PBH 
density is completely negligible (7]]2l[] ; the fraction of the 
energy density of the universe in PBHs at the time they 
form is given by p5(| 



Ppbh 

Ptot 



cr(AfH) exp 
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18ct 2 (M h ) 
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If n = 1 then a(M n ) = ct(M h ,o) = 9.5 x 1CT 5 and 
(3 ~ io(- 6xl ° 6 )! The Af~ 5 / 2 PBH mass function is there- 
fore cosmologically irrelevant. To produce an interesting 
density of PBHs either a power-law power-spectrum with 
more power on short scales (n > 1), or a spike in the den- 
sity perturbations spectrum is required @J|^] . In both of 
these cases if all PBHs which form at a given epoch have 
the same mass then the mass function is very sharply 
peaked. 

Niemeyer and Jedamzik jl6| have shown using numer- 
ical simulations that, as a consequence of near critical 
gravitational collapse Jij]] , at a fixed epoch PBHs with a 
range of masses form. The PBH mass is determined by 
the size of the fluctuation from which it formed: 



M BH = fcM H (<5 - 6 C )~ 
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P (6(M K ))dS(M K ) = 
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where <t(A/h) is the mass variance evaluated at horizon 
crossing [Q. For power law power spectra, P(k) oc fc™, 
where P(k) — (|<5k| 2 ) and n is the spectral index, 
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where 7, k and S c are constant for a given perturbation 
shape (for Mexican Hat shaped fluctuations 7 = 0.36, 
k = 2.85 and S c = 0.67.) 

It has been found that, for both power-law power spec- 
tra and flat spectra with a spike on a particular scale, 
in the limit where the number of PBHs formed is small 
enough to satisfy the observational constraints on their 
abundance, it can be assumed that all the PBHs form 
at a single horizon mass It is therefore possible to 
calculate the PBH initial number density per unit mass 
interval analytically 22 23]: 



where '0' and 'eq' denote quantities evaluated at 
the present day and matter-radiation equality, and 
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FIG. 1. The PBH initial number density, as a function of 
mass, for the Carr M BB 2 (solid line), Kim- Lee (dotted) and 
Niemeyer-Jedamzik (dashed) mass functions with parameters 
chosen such that the present day densities, ignoring evapora- 
tion, are the same (J7pbh = 1 x 10 -8 ). For clarity the Carr 
and Niemeyer-Jedamzik mass functions are multiplied by a 
factor of ten. 
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FIG. 2. The present day PBH number density, as a func- 
tion of mass, for the Carr M B ^ 2 (solid line), Kim-Lee (dot- 
ted) and Niemeyer-Jedamzik (dashed) mass functions with 
parameters chosen such that the present day number density 
of evaporating PBHs is 7V g = 4 x 10 3 pc -3 . 
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where ct(Mh) is defined in Eq. (g). We subsequently refer 
to this mass function as the Niemeyer-Jedamzik mass 
function. The physical number density of PBHs dilutes 
oc a~ 3 so that at any later time t 



d"BH , , _ dn-BH f T{t) 
dM BH [ ' ~ dAf BH V Ti 

where 7] can be related to Mm using 



(13) 



M m = M] 



HO 



T, 
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The three initial mass functions (Carr M Bn ' , Kim- 
Lee and Niemeyer-Jedamzik) are plotted in Fig. |l| with 
the parameters of each mass function chosen such that 
the energy density in PBHs is the same, corresponding to 
a present day PBH energy density of ^pbh = lx 10 -8 []. 
In Fig. ||we plot the present day mass functions with pa- 
rameters chosen such that J\f g — Ax 10 3 pc -3 . The Kim- 
Lee mass function, which arises from assuming that all 
PBHs which form at a given epoch have the same mass, 
is very sharply peaked. The Niemeyer-Jedamzik is far 
broader, with a long tail of low mass PBHs which would 
have evaporated since z ~ 700 and would contribute to 
the diffuse gamma-ray background. 



III. PBH EVAPORATION RATE 

The abundance of PBHs evaporating today is con- 
strained by their contribution to the diffuse gamma-ray 
background Jl^,[iJQ|l],|3) . All PBHs which have evap- 
orated since z = 700 contribute to the diffuse gamma-ray 
background |^5| , therefore the resulting constraint on the 
number of PBHs evaporating today depends on the PBH 
mass function; the broader the mass function the tighter 
the constraint on the number of PBHs evaporating to- 
day. The diffuse gamma-ray bound has been recalcu- 
lated using recent measurements by COMPTEL |2(J and 
EGRET [^7j for the Kim-Lee mass function p j j and also 
for the Niemeyer-Jedamzik mass function ]23|. 



Since the Kim-Lee and Carr M, 



-5/2 
BH 



mass function 



both have a sharp lower cut-off, the diffuse gamma- 
ray constraint on the present day PBH evaporation rate 
is the same for both mass functions [|] , with the re- 
cent COMPTEL and EGRET measurements tightening 
the constraints on Af g by a factor of about 1.5 [E8L 
For the Niemeyer-Jedamzik mass function the relevant 
constraint on PBHs with mass M — M* is <t(Mh = 
1.54 x 10 15 g) < 0.056. Inserting this in Eq. (|lg) leads 
to Afg < 2.7 x 10~ 4 /i 2 pc~ 3 so that 



dn B H a(M*) 



dt 



MS 



M w 5.2 x IQ- U h 2 r] pc~ 3 yr^ 1 . (15) 



This calculation assumes that the cosmological constant 
is zero, Qa = 0. If fl\ = 0.7, as current observation 



'Strictly speaking in the case of the Niemeyer-Jedamzik 
mass function the present day density will be less than this 
as the low mass tail of PBHs with Mbh < 5 x 10 14 g will have 
evaporated 

^The resulting constraint on the mass density in PBHs of all 
masses, Opbh, is far tighter for the Kim-Lee mass function 
however, since that mass function is far more sharply peaked. 
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appear to indicate [E9j , then the constraint is tightened to 
cr(M*) < 0.023 [|| and A" g < 4.5 x 10" 157 /i 2 pc~ 3 so that 
the present day rate of PBH evaporation is completely 
negligible. 

We caution that the constraints on Af g are very sensi- 
tive to small changes in the limits on <t(Mh) (or equiva- 
lently n) due to the exponential factor which arises in the 
expressions for dnBH/dMeH (Eqs. @ and (|l2|)). This 
illustrates the fact that to produce an interesting (i.e. 
non-negligible) density of PBHs (evaporating today or 
otherwise) requires fine tuning of the size of the den- 
sity perturbations from which the they form. We have 
also seen that if the PBH mass function has significant 
finite width then the diffuse gamma-ray bound (which 
constrains the abundance of all PBHs which have evapo- 
rated since a redshift of z = 700) places far tighter limits 
on the present day rate of PBH evaporation than if the 
PBH mass function is sharply peaked. 



IV. PBH CONCENTRATION 

The local density enhancement factor is usually cal- 
culated |H by assuming an isothermal halo so that the 
density at galactocentric radius R (outside the core ra- 
dius R c ) is related to the asymptotic circular velocity 
Voo ~ 220kms- 1 by M 
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V- 



AttGR 2 
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(16) 



which is in broad agreement with the value calculated 
from the asymptotic circular velocity assuming that the 
halo is isothermal. Since evaporating PBHs could only be 
detected by B ATSE out to distances of order a parsec || , 
it appears that, if the PBH mass function has significant 
width as recent numerical simulations suggest, then the 
local density enhancement factor is not large enough to 
produce a local PBH evaporation rate comparable with 
the observed frequency of short period gamma-ray bursts 
(10-20yr- 1 ). 



V. CONCLUSIONS 

Calculations of the present PBH evaporation rate have 
traditionally assumed that the PBH mass function varies 
- M B h /2 . This mass function only arises if the density 
perturbations from which the PBHs form have a scale 
invariant power spectrum, in which case the PBH density 
is completely negligible (7|j2l|]. We have recalculated the 
present PBH evaporation rate, using the bounds which 
arise from the diffuse gamma-ray background |^,^3| , for 
the sharply peaked mass function which arises if all PBHs 
which form at a given epoch have the same mass and 
also for the broader PBH mass function found by recent 
numerical simulations. We find that if the PBH mass 
function has significant finite width it is not possible to 
produce a present day PBH evaporation rate comparable 
with the observed short period gamma-ray burst rate. 
This could also have implications for other attempts to 
detect evaporating PBHs. 



Taking our galactocentric radius to be Rq ■ 
8.5(±l.l)kpc gives a local density enhancement factor 



// * 4.r,( ±o.(i) \ :io 7/--' ( — 



(17) 



where J7h is the fraction of the critical density in galactic 
halos. Halzen et. al. |p| claim, without explicit calcula- 
tion, a larger value for the local density enhancement: 
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(18) 



which apparently |3l|] relies on the assumption that PBHs 
are concentrated to the same extent as luminous matter. 

An alternative procedure is to calculate r\ directly from 
estimates of the local halo density. The local halo density 
is poorly known, with the most recent estimates giving 
pi = 9-2±|;f x 10~ 25 gcm" 3 Q. Taking the currently 
favoured value of the present day mass density f2 m = 
0.3 [fH gives p m « 5.6 x lO" 30 ^ 2 gem" 3 . Combining 
these values produces 



Tj = 1.6_ 



0.8 



x 10 5 h~ 



03 



(19) 
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